Charge ordering and long-range interactions in layered transition metal oxides 
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We study the competition between long-range and short-range interactions among holes within the 
spin density wave picture of layered transition metal oxides. We focus on the problem of charge 
ordering and the charge phase diagram. We show that the main interactions are the long-range 
Coulomb interaction and a dipolar short-range interaction generated by the short-range antifer- 
romagnetic fluctuations. We flnd four different phases depending on the strength of the dipolar 
interaction and the density of holes: Wigner crystal, diagonal stripes, horizontal-vertical stripes 
(loops) and a glassy-clumped phase. We discuss the effect of temperature, disorder and lattice 
effects on these phases. 

PACS numbers: 71.10.Hf, 71.20.Be, 71.10.-w 



Recently there has been much interest in the charge or- 
dering and domain wall formation at mesoscopic scales in 
doped transition metal oxides Q . A popular example of 
such orderings are stripes, i.e., linear arrays of holes sep- 
arated by an antiferromagnetically (AF) ordered back- 
ground. The formation of domain walls and stripes has 
been discussed in terms of the proximity to phase separa- 
tion . Macroscopic phase separation has been observed 
in La2Cu04+5 and stripes have been observed exper- 
imentally in La2_KSrj:Ni04+y in many different exper- 
iments including direct high-resolution electron diffrac- 
tion Magnetic susceptibility measurements nu- 
clear quadrupole resonance |^ and muon spin resonance 
[0 indicate formation of domains in La2-xSra;Cu04. 
Stripes have also been seen in Lai-ajCa^jMnOa for specific 
commensurate values of doping A direct evidence 
for stripe formation was observed in neutron scattering 
in Lai.g_a;Ndo.4Sra;Cu04 |^]. Moreover, recent neutron 
scattering experiments are not inconsistent with stripe 
phases in other high temperature superconductors such 
as YBa2Cu307-5 0,0. 

On the theoretical side, stripes have been proposed as 
a result of a competition between short-range attractive 
interaction between holes from the breaking of AF bonds 
and the long-range Coulomb interaction Indirect 
support for this picture has been given in terms of the 
mapping of the problem into effective spin models [p^ . 
Striped phases have been obtained within mean field ap- 
proaches to the short-range Hubbard or t — J models 
which are only able to generate insulating states p4[ . 
Numerical methods in these models have not been able 
to confirm this picture except for recent Density Matrix 
Renormalization Group simulations (DMRG) Jisf . 

In this paper we present a numerical approach to the 
problem of holes moving in an AF insulator in the pres- 
ence of long-range Coulomb forces. The ability to han- 
dle long-range Coulomb interactions at finite density has 
been enhanced recently in the area of molecular physics: 
assuming a computational cell of arbitrary geometry and 



cyclic boundary conditions it is possible to sum inter- 
actions of particles with all of their images residing in 
cells obtained by translation from the original computa- 
tional cell |]l6| , |l7| . On making integral transformations. 
Coulomb interactions are computed by summing over 
fast-convergent Bessel functions with great accuracy. 

Using Monte Carlo (MC) and molecular dynamics 
(MD) methods, we systematically study the interplay be- 
tween long-range Coulomb interaction and short-range 
AF interactions of dipolar nature which we take to have 
both isotropic and anisotropic components (depending on 
the lattice structure). Our main result is summarized in 
the phase diagram of Figj^. In the absence of disorder we 
find four phases depending on the density of holes and 
the characteristic AF energy scales: a Wigner crystal, 
diagonal stripes, horizontal-vertical stripes (loops) and a 
glassy-clumped phase. The order parameter for charge 
ordering is the Fourier transform of the hole density: 



p(q) 
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where is the position of the i*'' hole and N is the total 
number of holes. A peak in p(q) at some wave-vector 
q = K indicates ordering. 

Our starting point is the spin density wave (SDW) pic- 
ture of the layered transition metal oxides which has been 
very successful in describing the insulating AF phase of 
these systems ||lj,|l^. In this picture the electrons move 
with hopping energy t in the self-consistent staggered 
field of its spin. Because the translational symmetry of 
the system is broken, the electronic band is split into 
upper and lower Hubbard bands pO| . These are sep- 
arated by the Mott-Hubbard gap. A, and at half fill- 
ing the lower band is filled and the upper one is empty. 
This picture is consistent with the angle resolved pho- 
toemission data in the layered AF insulator Sr2Cu02Cl2 
pl[ . By doping the system with holes with planar den- 
sity Os and at low temperatures, T {ksT << A), we 
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focus entirely on the lower band which has a maximum 
at k = Q/2 = (±l,±l)7r/(2a), where a is the lattice 
spacing. It can be shown that the holes interact via two 
different mechanisms: a short-range attractive force due 
to AF bond breaking and a long-range dipolar interac- 
tion due to the distortion of the AF background | p^ , p2[ . 
It was shown that this dipolar interaction gives rise to 
spiral distortions of the AF background [^,^. The 
dipole moment associated with each hole is due to the 
virtual hopping of holes between neighboring sites and 
scales with the AF magnetic energy. The dipolar inter- 
action between two holes with dipole moments di 2 at 
distance r apart has the form: 



(di.d2)-^(di.r)(d2T) 



(2) 



which is rotationally invariant. It is also possible to show 
using Ward identities that the spin part of the problem 
can be described by a two dimensional (2D) non-linear a 
model in the long wavelength limit (24j. At finite T the 
system is magnetically disordered and characterized by a 
finite magnetic correlation length, ^ |25|. Thus, at finite 
T the dipolar interaction between the holes, mediated by 
the antiferromagnet, is actually short-ranged. However, 
besides the AF interactions the holes also feel the long- 
range Coulomb interaction. This is clear if we consider 
that rs = ro/ao (where ro is the mean inter-particle dis- 
tance and oq is the Bohr radius) is very large in the under- 
doped systems (r^ « 8). Thus, the interaction energy be- 
tween the holes, which behaves like e^/(ao?'s), is certainly 
more important than the kinetic term (« e^/(aor^)) at 
low densities. This implies that the interaction terms 
should be treated first and the kinetic energy as a pertur- 
bation. Finally, each hole carries a spin degree of freedom 
as well, but it is possible to show that the overall spin 
energy is minimized in the spin anti-symmetric channel, 
as we assume here. Thus, in our approach, we are left 
with only the charge channel and the interaction between 
two holes, 1 and 2, has the form (see Eq. (||)) 
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where q is the hole charge, 9 is the angle made between r 
and a fixed axis and 4>i^2 are the angles of the dipoles rela- 
tive to the same fixed axis. A is the strength of the short- 
range bond-breaking interaction and B is the strength of 
the dipolar interaction, which we will assume to be inde- 
pendent variables. The magnetic correlation length ^ is 
obtained from neutron scattering measurements pGf . It 
is also worth mentioning that the 1/2-filled Landau level 
problem has been mapped into an interacting 2D dipole 
gas 1^ . More recently it was shown that the same type 
of description is possible for a 2D electron gas even in the 
absence of a magnetic field . 

In general, the many-body problem of holes in an AF 
background is extremely complicated, involving many- 



particle interaction terms. However, at low densities it is 
reasonable to assume that the interaction of any two holes 
is weakly perturbed by other holes, and the total poten- 
tial energy can be expressed in terms of two-particle en- 
ergies. Therefore, in our numerical calculations we study 
the physics of N holes interacting via V{r) as given in 
(^). We assume a rectangular computational box of size 
Lx X Ly with Lx, Ly up to 100 unit cells in a Cu02 plane. 
At the beginning of each simulation we place the holes 
at random and assign to each hole a dipole moment of 
constant size, but random direction. We find a minimum 
of the total potential in this system using three different 
methods: MC method, Langevin MD and a hybrid MC- 
MD method . All three methods yield essentially the 
same results. Since the system exhibits several phases 
(see Fig. ^ for some values of the input parameters, its 
ground state is not always well defined and may, in fact, 
depend on the initial and boundary conditions. Hence, 
in order to rapidly reach a hole configuration with the 
lowest global minimum energy we perform simulated an- 
nealing from high temperatures. 

For _B = we find the Wigner crystal with small dis- 
tortions to be the state of lowest energy, as expected . 
The small distortion of the crystal structure is due to the 
periodicity, which introduces a small spatial anisotropy 
into the system due to the rectangular shape of the com- 
putational box. Increasing A while retaining i? = re- 
duces the lattice constant of the Wigner crystal until a 
critical value is reached where holes group together. For 
A = and finite B the situation is quite different. At 
small B and larger densities the Wigner crystal is unsta- 
ble and a new phase with diagonal stripes is formed. This 
phase is characterized by ferro-dipolar order (see Fig. ^). 
The situation here is very similar to that observed in 
La2-2;Sra;Ni04+y As shown in Fig. at larger val- 
ues of -B a line stripe is formed, which, with increasing 
density tends to close into loops, forming a checkerboard 
pattern. Importantly, the loop formation is accompanied 
by dipole orientation along the straight portion of a loop 
with gradual rotation by 7r/2 at each corner Due to 
the rotation of dipoles at corners the loops interact, and 
eventually form the checkerboard pattern The size 
of the inter-hole distance within a line is determined by 
the ratio of B and the Coulomb energy; the loop sizes 
are determined by the hole density alone. These results 
appear to be consistent with the DMRG solution of the 
t-J model ||l5[. If B is increased further the dipolar in- 
teraction becomes dominant over the average Coulomb 
interaction; the well-defined pattern disappears and one 
observes star shaped clumps of holes, which can, at suf- 
ficiently high density, form another geometric structure 
(e.g., a Wigner crystal of clumps). We remark that in 
all phases a non- vanishing value of A leads to a decrease 
in the effective value of B at which the transitions oc- 
cur (Fig. ^; the isotropic term A alone never produces 
any non-trivial geometric phase (e.g., stripes), even with 
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inclusion of lattice effects. We find that the transition be- 
tween the ferro-dipolar and the stripe phase is first order, 
while other transitions appear to be of second order . 
The stripe tension of the hole patterns will be quantified 
elsewhere Q. 

In the cases presented above we have assumed uniform 
dipolar interaction. It is well known that there are or- 
thorhombic and tetragonal distortions in practically all 
transition metal oxides. In particular static stripe for- 
mation has only been observed in the low temperature 
tetragonal phase of Lai.6-a;Ndo.4Sra;Cu04 Q. In order 
to study the influence of the anisotropy we assume that 
the dipole sizes along x and y directions have anisotropy 
a {a = 1 corresponds to the isotropic case). Figure ||c 
shows our solution for a = 0.8: the symmetry is bro- 
ken and a stripe superlattice is formed, with a charge 
ordering vector K = (7r/^)x, where £ is the inter-stripe 
distance. In the SDW model the Fourier transform of the 
magnetization SYq) = (5*2 (q)) is slaved to such that 
a peak at K in (|l|) leads to a peak at Q ± K in S{q) Q . 
Thus our results yield a neutron peak at (7r/a±7r/^, tt/o). 
Assuming twinning, this would imply neutron peaks at 
(jr/aztTr/i, 7r/a±7r/^) in agreement with experiment 
The same is obtained in the checkerboard phase (see Figs. 
|cand|d). If one includes the kinetic energy |p3| , instead 
of static stripe formation one would obtain dynamical 
stripes like those believed to exist in La2Sra;Cu04. In 
this case the Fermi surface of the system should be mod- 
ified by the superlattice formation p^ . 

Our results are somewhat sensitive to the applied 
boundary conditions: first, the exact size of the checker- 
boards depends on its commensuration with the com- 
putational box, which, in turn depends on the density. 
On increasing of the size of the cornputational box, the 
checkerboard pattern shown in Fig. acquires point or 
line defects . This leads to the reduction in the higher 
order peaks observed in Fig. |^d with no change in their 
wave-numbers. Second, in a finite system with appropri- 
ate charge background, the holes do not form geometric 
phases, although they still form stripes js^. However, in 
this case even a very small anisotropy (a ^ 0.95) again 
leads to stripe formation as in Fig. ^ [P3| . 

We have also studied impurity effects (from defects or 
charged counter- ions). For example, we place the same 
number of impurities as holes, randomly in a plane a dis- 
tance d = 6 A above the plane to simulate the situation in, 
e.g., Sr doped cuprates and consider the unscreened at- 
tractive Coulomb interaction between impurity and hole. 
The charge pattern produced is very sensitive to impurity 
doping (see Fig. |2|f). The Wigner crystal becomes glassy 
with no obvious sign of charge ordering. This happens 
because the attractive Coulomb energy between impuri- 
ties and holes scales like e^/rf while the average inter- 
hole Coulomb energy behaves like ^fWl. Thus when 
<Js < the holes are pinned by impurities. Most 

strikingly, all other phases are unstable towards finite 



stripe formation. The loops and diagonal stripes tend 
to deform to pass very close to the impurities in order 
to maximize the attractive energy. However, the dipole 
interaction is sufficient to retain the main orientation. 
This leads us to conjecture that with the addition of the 
kinetic energy the holes can move in string segments in 
an orientation given basically by the phase diagram of 
the clean system. These string segments are kept to- 
gether by the dipolar interaction (i.e., string tension). 
The stripe motion would then be caused by mesoscopic 
thermal or quantum tunneling of the finite strings be- 
tween the minima of the overall potential. This would 
lead to non-linear effects in the low temperature field de- 
pendent conductivity |^ and unusual T dependence of 
the conductivity |3^. We have also performed simula- 
tions in the presence of a realistic underlying periodic 
lattice and have found that this creates slight distortions 
in the phases, pinning loops more strongly |Q . In partic- 
ular, the peaks in p(q) sharpen in some of these phases. 
Finally, at finite T melting of the phases occurs because 
of the small energy scales and large entropy in long-range 
Coulomb tails. 

In summary, using a novel numerical technique, we 
have studied the competition between long-range and 
short-range interactions and its impact on hole order- 
ing in layered transition metal oxides. Employing the 
SDW picture of these systems, we have studied the short- 
range attractive force and the dipolar force generated 
by the short-range AF fluctuations together with long- 
range Coulomb forces for a 2D layer. We have found a 
rich phase diagram for the clean system which includes 
a Wigner solid, stripes, loops and a glassy phase. This 
phase diagram is consistent with several different exper- 
imental measurements. We have also found this system 
to be rather sensitive to the presence of charged impuri- 
ties. However, the stripe phases survive as finite stripe 
segments which we believe is key to understanding these 
systems. Finally we have also investigated the effects of 
a periodic lattice, which leads to minor distortions of the 
various phases. These results are of great importance 
for the understanding of the phase diagram of layered 
transition metal oxides. 
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FIG. 1. Phase diagram as a function of the hole density 
and the strength of the dipolar interaction, B, for A — 0. 
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FIG. 2. Geometric phases resulting from the competition 
of dipolar and Coulomb interactions: panels (a) and (c) show 
holes (open circles) with their dipole orientations, in a small 
section of the computational box, for the ferro-dipolar and 
stripe phases, respectively; panels (b) and (d) show contour 
plots of the hole density in momentum space (see Eq. (1)) for 
the phases shown in panels (a) and (c). Panel (e) shows the 
stripe phase obtained with dipole anisotropy of 0.8, and panel 
(f) hole positions in the presence of impurities (crosses). 
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